Reprogramming of the developing heart by Hif1a-deficient sympathetic system and maternal diabetes exposure.

Introduction: Maternal diabetes is a recognized risk factor for both short-term and long-term complications in offspring. Beyond the direct teratogenicity of maternal diabetes, the intrauterine environment can influence the offspring's cardiovascular health. Abnormalities in the cardiac sympathetic system are implicated in conditions such as sudden infant death syndrome, cardiac arrhythmic death, heart failure, and certain congenital heart defects in children from diabetic pregnancies. However, the mechanisms by which maternal diabetes affects the development of the cardiac sympathetic system and, consequently, heightens health risks and predisposes to cardiovascular disease remain poorly understood. Methods and results: In the mouse model, we performed a comprehensive analysis of the combined impact of a Hif1a-deficient sympathetic system and the maternal diabetes environment on both heart development and the formation of the cardiac sympathetic system. The synergic negative effect of exposure to maternal diabetes and Hif1a deficiency resulted in the most pronounced deficit in cardiac sympathetic innervation and the development of the adrenal medulla. Abnormalities in the cardiac sympathetic system were accompanied by a smaller heart, reduced ventricular wall thickness, and dilated subepicardial veins and coronary arteries in the myocardium, along with anomalies in the branching and connections of the main coronary arteries. Transcriptional profiling by RNA sequencing (RNA-seq) revealed significant transcriptome changes in Hif1a-deficient sympathetic neurons, primarily associated with cell cycle regulation, proliferation, and mitosis, explaining the shrinkage of the sympathetic neuron population. Discussion: Our data demonstrate that a failure to adequately activate the HIF-1α regulatory pathway, particularly in the context of maternal diabetes, may contribute to abnormalities in the cardiac sympathetic system. In conclusion, our findings indicate that the interplay between deficiencies in the cardiac sympathetic system and subtle structural alternations in the vasculature, microvasculature, and myocardium during heart development not only increases the risk of cardiovascular disease but also diminishes the adaptability to the stress associated with the transition to extrauterine life, thus increasing the risk of neonatal death.

Human bone marrow stromal cells: the impact of anticoagulants on stem cell properties.

Background: Bone marrow stromal cells (BMSCs) are the source of multipotent stem cells, which are important for regenerative medicine and diagnostic purposes. The isolation of human BMSCs from the bone marrow (BM) cavity using BM aspiration applies the method with collection into tubes containing anticoagulants. Interactions with anticoagulants may affect the characteristics and composition of isolated BMSCs in the culture. Thus, we investigated how anticoagulants in isolation procedures and cultivation affect BMSC molecular characteristics. Methods: BM donors (age: 48-85 years) were recruited from the hematology clinic. BM aspirates were obtained from the iliac crest and divided into tubes coated with ethylenediaminetetraacetic acid (EDTA) or heparin anticoagulants. Isolated BMSCs were analyzed by flow cytometry and RNA-seq analysis. Further cellular and molecular characterizations of BMSCs including CFU, proliferation and differentiation assays, cytometry, bioenergetic assays, metabolomics, immunostaining, and RT-qPCR were performed. Results: The paired samples of isolated BMSCs obtained from the same patient showed increased cellular yield in heparin vs. EDTA samples, accompanied by the increased number of CFU colonies. However, no significant changes in molecular characteristics were found between heparin- and EDTA-isolated BMSCs. On the other hand, RNA-seq analysis revealed an increased expression of genes involved in nucleotide metabolism and cellular metabolism in cultivated vs. non-cultivated BMSCs regardless of the anticoagulant, while genes involved in inflammation and chromatin remodeling were decreased in cultivated vs. non-cultivated BMSCs. Conclusion: The type of anticoagulant in BMSC isolation did not have a significant impact on molecular characteristics and cellular composition, while in vitro cultivation caused the major change in the transcriptomics of BMSCs, which is important for future protocols using BMSCs in regenerative medicine and clinics.

NEUROD1 reinforces endocrine cell fate acquisition in pancreatic development.

NEUROD1 is a transcription factor that helps maintain a mature phenotype of pancreatic ß cells. Disruption of Neurod1 during pancreatic development causes severe neonatal diabetes; however, the exact role of NEUROD1 in the differentiation programs of endocrine cells is unknown. Here, we report a crucial role of the NEUROD1 regulatory network in endocrine lineage commitment and differentiation. Mechanistically, transcriptome and chromatin landscape analyses demonstrate that Neurod1 inactivation triggers a downregulation of endocrine differentiation transcription factors and upregulation of non-endocrine genes within the Neurod1-deficient endocrine cell population, disturbing endocrine identity acquisition. Neurod1 deficiency altered the H3K27me3 histone modification pattern in promoter regions of differentially expressed genes, which resulted in gene regulatory network changes in the differentiation pathway of endocrine cells, compromising endocrine cell potential, differentiation, and functional properties.

Dysregulation of hypoxia-inducible factor 1α in the sympathetic nervous system accelerates diabetic cardiomyopathy.

BACKGROUND: An altered sympathetic nervous system is implicated in many cardiac pathologies, ranging from sudden infant death syndrome to common diseases of adulthood such as hypertension, myocardial ischemia, cardiac arrhythmias, myocardial infarction, and heart failure. Although the mechanisms responsible for disruption of this well-organized system are the subject of intensive investigations, the exact processes controlling the cardiac sympathetic nervous system are still not fully understood. A conditional knockout of the Hif1a gene was reported to affect the development of sympathetic ganglia and sympathetic innervation of the heart. This study characterized how the combination of HIF-1α deficiency and streptozotocin (STZ)-induced diabetes affects the cardiac sympathetic nervous system and heart function of adult animals. METHODS: Molecular characteristics of Hif1a deficient sympathetic neurons were identified by RNA sequencing. Diabetes was induced in Hif1a knockout and control mice by low doses of STZ treatment. Heart function was assessed by echocardiography. Mechanisms involved in adverse structural remodeling of the myocardium, i.e. advanced glycation end products, fibrosis, cell death, and inflammation, was assessed by immunohistological analyses. RESULTS: We demonstrated that the deletion of Hif1a alters the transcriptome of sympathetic neurons, and that diabetic mice with the Hif1a-deficient sympathetic system have significant systolic dysfunction, worsened cardiac sympathetic innervation, and structural remodeling of the myocardium. CONCLUSIONS: We provide evidence that the combination of diabetes and the Hif1a deficient sympathetic nervous system results in compromised cardiac performance and accelerated adverse myocardial remodeling, associated with the progression of diabetic cardiomyopathy.

Cortical glia in SOD1(G93A) mice are subtly affected by ALS-like pathology.

The role of glia in amyotrophic lateral sclerosis (ALS) is undeniable. Their disease-related activity has been extensively studied in the spinal cord, but only partly in the brain. We present herein a comprehensive study of glia in the cortex of SOD1(G93A) mice-a widely used model of ALS. Using single-cell RNA sequencing (scRNA-seq) and immunohistochemistry, we inspected astrocytes, microglia, and oligodendrocytes, in four stages of the disease, respecting the factor of sex. We report minimal changes of glia throughout the disease progression and regardless of sex. Pseudobulk and single-cell analyses revealed subtle disease-related transcriptional alterations at the end-stage in microglia and oligodendrocytes, which were supported by immunohistochemistry. Therefore, our data support the hypothesis that the SOD1(G93A) mouse cortex does not recapitulate the disease in patients, and we recommend the use of a different model for future studies of the cortical ALS pathology.

Complement C3a treatment accelerates recovery after stroke via modulation of astrocyte reactivity and cortical connectivity.

Despite advances in acute care, ischemic stroke remains a major cause of long-term disability. Approaches targeting both neuronal and glial responses are needed to enhance recovery and improve long-term outcome. The complement C3a receptor (C3aR) is a regulator of inflammation with roles in neurodevelopment, neural plasticity, and neurodegeneration. Using mice lacking C3aR (C3aR-/-) and mice overexpressing C3a in the brain, we uncovered 2 opposing effects of C3aR signaling on functional recovery after ischemic stroke: inhibition in the acute phase and facilitation in the later phase. Peri-infarct astrocyte reactivity was increased and density of microglia reduced in C3aR-/- mice; C3a overexpression led to the opposite effects. Pharmacological treatment of wild-type mice with intranasal C3a starting 7 days after stroke accelerated recovery of motor function and attenuated astrocyte reactivity without enhancing microgliosis. C3a treatment stimulated global white matter reorganization, increased peri-infarct structural connectivity, and upregulated Igf1 and Thbs4 in the peri-infarct cortex. Thus, C3a treatment from day 7 after stroke exerts positive effects on astrocytes and neuronal connectivity while avoiding the deleterious consequences of C3aR signaling during the acute phase. Intranasal administration of C3aR agonists within a convenient time window holds translational promise to improve outcome after ischemic stroke.

Glioblastoma and cerebral organoids: development and analysis of an in vitro model for glioblastoma migration.

It is currently challenging to adequately model the growth and migration of glioblastoma using two-dimensional (2D) in vitro culture systems as they quickly lose the original, patient-specific identity and heterogeneity. However, with the advent of three-dimensional (3D) cell cultures and human-induced pluripotent stem cell (iPSC)-derived cerebral organoids (COs), studies demonstrate that the glioblastoma-CO (GLICO) coculture model helps to preserve the phenotype of the patient-specific tissue. Here, we aimed to set up such a model using mature COs and develop a pipeline for subsequent analysis of cocultured glioblastoma. Our data demonstrate that the growth and migration of the glioblastoma cell line within the mature COs are significantly increased in the presence of extracellular matrix proteins, shortening the time needed for glioblastoma to initiate migration. We also describe in detail the method for the visualization and quantification of these migrating cells within the GLICO model. Lastly, we show that this coculture model (and the human brain-like microenvironment) can significantly transform the gene expression profile of the established U87 glioblastoma cell line into proneural and classical glioblastoma cell types.

Recent advances in deciphering oligodendrocyte heterogeneity with single-cell transcriptomics.

Oligodendrocytes (OL) have been for decades considered a passive, homogenous population of cells that provide support to neurons, and show a limited response to pathological stimuli. This view has been dramatically changed by the introduction of powerful transcriptomic methods that have uncovered a broad spectrum of OL populations that co-exist within the healthy central nervous system (CNS) and also across a variety of diseases. Specifically, single-cell and single-nucleus RNA-sequencing (scRNA-seq, snRNA-seq) have been used to reveal OL variations in maturation, myelination and immune status. The newly discovered immunomodulatory role suggests that OL may serve as targets for future therapies. In this review, we summarize the current understanding of OL heterogeneity in mammalian CNS as revealed by scRNA-seq and snRNA-seq. We provide a list of key studies that identify consensus marker genes defining the currently known OL populations. This resource can be used to standardize analysis of OL related datasets and improve their interpretation, ultimately leading to a better understanding of OL functions in health and disease.

Evolutionary conservation of maternal RNA localization in fishes and amphibians revealed by TOMO-Seq.

Asymmetrical localization of biomolecules inside the egg, results in uneven cell division and establishment of many biological processes, cell types and the body plan. However, our knowledge about evolutionary conservation of localized transcripts is still limited to a few models. Our goal was to compare localization profiles along the animal-vegetal axis of mature eggs from four vertebrate models, two amphibians (Xenopus laevis, Ambystoma mexicanum) and two fishes (Acipenser ruthenus, Danio rerio) using the spatial expression method called TOMO-Seq. We revealed that RNAs of many known important transcripts such as germ layer determinants, germ plasm factors and members of key signalling pathways, are localized in completely different profiles among the models. It was also observed that there was a poor correlation between the vegetally localized transcripts but a relatively good correlation between the animally localized transcripts. These findings indicate that the regulation of embryonic development within the animal kingdom is highly diverse and cannot be deduced based on a single model.

Mitochondrial respiration supports autophagy to provide stress resistance during quiescence.

Mitochondrial oxidative phosphorylation (OXPHOS) generates ATP, but OXPHOS also supports biosynthesis during proliferation. In contrast, the role of OXPHOS during quiescence, beyond ATP production, is not well understood. Using mouse models of inducible OXPHOS deficiency in all cell types or specifically in the vascular endothelium that negligibly relies on OXPHOS-derived ATP, we show that selectively during quiescence OXPHOS provides oxidative stress resistance by supporting macroautophagy/autophagy. Mechanistically, OXPHOS constitutively generates low levels of endogenous ROS that induce autophagy via attenuation of ATG4B activity, which provides protection from ROS insult. Physiologically, the OXPHOS-autophagy system (i) protects healthy tissue from toxicity of ROS-based anticancer therapy, and (ii) provides ROS resistance in the endothelium, ameliorating systemic LPS-induced inflammation as well as inflammatory bowel disease. Hence, cells acquired mitochondria during evolution to profit from oxidative metabolism, but also built in an autophagy-based ROS-induced protective mechanism to guard against oxidative stress associated with OXPHOS function during quiescence.Abbreviations: AMPK: AMP-activated protein kinase; AOX: alternative oxidase; Baf A: bafilomycin A1; CI, respiratory complexes I; DCF-DA: 2',7'-dichlordihydrofluorescein diacetate; DHE: dihydroethidium; DSS: dextran sodium sulfate; ΔΨmi: mitochondrial inner membrane potential; EdU: 5-ethynyl-2'-deoxyuridine; ETC: electron transport chain; FA: formaldehyde; HUVEC; human umbilical cord endothelial cells; IBD: inflammatory bowel disease; LC3B: microtubule associated protein 1 light chain 3 beta; LPS: lipopolysaccharide; MEFs: mouse embryonic fibroblasts; MTORC1: mechanistic target of rapamycin kinase complex 1; mtDNA: mitochondrial DNA; NAC: N-acetyl cysteine; OXPHOS: oxidative phosphorylation; PCs: proliferating cells; PE: phosphatidylethanolamine; PEITC: phenethyl isothiocyanate; QCs: quiescent cells; ROS: reactive oxygen species; PLA2: phospholipase A2, WB: western blot.

Isolation and Characterization of Highly Pure Type A Spermatogonia From Sterlet (Acipenser ruthenus) Using Flow-Cytometric Cell Sorting.

Sturgeons are among the most ancient linages of actinopterygians. At present, many sturgeon species are critically endangered. Surrogate production could be used as an affordable and a time-efficient method for endangered sturgeons. Our study established a method for identifying and isolating type A spermatogonia from different developmental stages of testes using flow cytometric cell sorting (FCM). Flow cytometric analysis of a whole testicular cell suspension showed several well-distinguished cell populations formed according to different values of light scatter parameters. FCM of these different cell populations was performed directly on glass slides for further immunocytochemistry to identify germ cells. Results showed that the cell population in gate P1 on a flow cytometry plot (with high forward scatter and high side scatter parameter values) contains the highest amount of type A spermatogonia. The sorted cell populations were characterized by expression profiles of 10 germ cell specific genes. The result confirmed that setting up for the P1 gate could precisely sort type A spermatogonia in all tested testicular developmental stages. The P2 gate, which was with lower forward scatter and side scatter values mostly, contained type B spermatogonia at a later maturing stage. Moreover, expressions of plzf, dnd, boule, and kitr were significantly higher in type A spermatogonia than in later developed germ cells. In addition, plzf was firstly found as a reliable marker to identify type A spermatogonia, which filled the gap of identification of spermatogonial stem cells in sterlet. It is expected to increase the efficiency of germ stem cell culture and transplantation with plzf identification. Our study thus first addressed a phenotypic characterization of a pure type A spermatogonia population in sterlet. FCM strategy can improve the production of sturgeons with surrogate broodstock and further the analysis of the cellular and molecular mechanisms of sturgeon germ cell development.

Transient astrocyte-like NG2 glia subpopulation emerges solely following permanent brain ischemia.

NG2 glia display wide proliferation and differentiation potential under physiological and pathological conditions. Here, we examined these two features following different types of brain disorders such as focal cerebral ischemia (FCI), cortical stab wound (SW), and demyelination (DEMY) in 3-month-old mice, in which NG2 glia are labeled by tdTomato under the Cspg4 promoter. To compare NG2 glia expression profiles following different CNS injuries, we employed single-cell RT-qPCR and self-organizing Kohonen map analysis of tdTomato-positive cells isolated from the uninjured cortex/corpus callosum and those after specific injury. Such approach enabled us to distinguish two main cell populations (NG2 glia, oligodendrocytes), each of them comprising four distinct subpopulations. The gene expression profiling revealed that a subpopulation of NG2 glia expressing GFAP, a marker of reactive astrocytes, is only present transiently after FCI. However, following less severe injuries, namely the SW and DEMY, subpopulations mirroring different stages of oligodendrocyte maturation markedly prevail. Such injury-dependent incidence of distinct subpopulations was also confirmed by immunohistochemistry. To characterize this unique subpopulation of transient astrocyte-like NG2 glia, we used single-cell RNA-sequencing analysis and to disclose their basic membrane properties, the patch-clamp technique was employed. Overall, we have proved that astrocyte-like NG2 glia are a specific subpopulation of NG2 glia emerging transiently only following FCI. These cells, located in the postischemic glial scar, are active in the cell cycle and display a current pattern similar to that identified in cortical astrocytes. Astrocyte-like NG2 glia may represent important players in glial scar formation and repair processes, following ischemia.

Restriction of an intron size en route to endothermy.

Ca2+-insensitive and -sensitive E1 subunits of the 2-oxoglutarate dehydrogenase complex (OGDHC) regulate tissue-specific NADH and ATP supply by mutually exclusive OGDH exons 4a and 4b. Here we show that their splicing is enforced by distant lariat branch points (dBPs) located near the 5' splice site of the intervening intron. dBPs restrict the intron length and prevent transposon insertions, which can introduce or eliminate dBP competitors. The size restriction was imposed by a single dominant dBP in anamniotes that expanded into a conserved constellation of four dBP adenines in amniotes. The amniote clusters exhibit taxon-specific usage of individual dBPs, reflecting accessibility of their extended motifs within a stable RNA hairpin rather than U2 snRNA:dBP base-pairing. The dBP expansion took place in early terrestrial species and was followed by a uridine enrichment of large downstream polypyrimidine tracts in mammals. The dBP-protected megatracts permit reciprocal regulation of exon 4a and 4b by uridine-binding proteins, including TIA-1/TIAR and PUF60, which promote U1 and U2 snRNP recruitment to the 5' splice site and BP, respectively, but do not significantly alter the relative dBP usage. We further show that codons for residues critically contributing to protein binding sites for Ca2+ and other divalent metals confer the exon inclusion order that mirrors the Irving-Williams affinity series, linking the evolution of auxiliary splicing motifs in exons to metallome constraints. Finally, we hypothesize that the dBP-driven selection for Ca2+-dependent ATP provision by E1 facilitated evolution of endothermy by optimizing the aerobic scope in target tissues.

NormQ: RNASeq normalization based on RT-qPCR derived size factors.

The merit of RNASeq data relies heavily on correct normalization. However, most methods assume that the majority of transcripts show no differential expression between conditions. This assumption may not always be correct, especially when one condition results in overexpression. We present a new method (NormQ) to normalize the RNASeq library size, using the relative proportion observed from RT-qPCR of selected marker genes. The method was compared against the popular median-of-ratios method, using simulated and real-datasets. NormQ produced more matches to differentially expressed genes in the simulated dataset and more distribution profile matches for both simulated and real datasets.

Decoding the Transcriptional Response to Ischemic Stroke in Young and Aged Mouse Brain.

Ischemic stroke is a well-recognized disease of aging, yet it is unclear how the age-dependent vulnerability occurs and what are the underlying mechanisms. To address these issues, we perform a comprehensive RNA-seq analysis of aging, ischemic stroke, and their interaction in 3- and 18-month-old mice. We assess differential gene expression across injury status and age, estimate cell type proportion changes, assay the results against a range of transcriptional signatures from the literature, and perform unsupervised co-expression analysis, identifying modules of genes with varying response to injury. We uncover downregulation of axonal and synaptic maintenance genetic program, and increased activation of type I interferon (IFN-I) signaling following stroke in aged mice. Together, these results paint a picture of ischemic stroke as a complex age-related disease and provide insights into interaction of aging and stroke on cellular and molecular level.

The role of nitric oxide during embryonic wound healing.

BACKGROUND: The study of the mechanisms controlling wound healing is an attractive area within the field of biology, with it having a potentially significant impact on the health sector given the current medical burden associated with healing in the elderly population. Healing is a complex process and includes many steps that are regulated by coding and noncoding RNAs, proteins and other molecules. Nitric oxide (NO) is one of these small molecule regulators and its function has already been associated with inflammation and angiogenesis during adult healing. RESULTS: Our results showed that NO is also an essential component during embryonic scarless healing and acts via a previously unknown mechanism. NO is mainly produced during the early phase of healing and it is crucial for the expression of genes associated with healing. However, we also observed a late phase of healing, which occurs for several hours after wound closure and takes place under the epidermis and includes tissue remodelling that is dependent on NO. We also found that the NO is associated with multiple cellular metabolic pathways, in particularly the glucose metabolism pathway. This is particular noteworthy as the use of NO donors have already been found to be beneficial for the treatment of chronic healing defects (including those associated with diabetes) and it is possible that its mechanism of action follows those observed during embryonic wound healing. CONCLUSIONS: Our study describes a new role of NO during healing, which may potentially translate to improved therapeutic treatments, especially for individual suffering with problematic healing.

Alternative assembly of respiratory complex II connects energy stress to metabolic checkpoints.

Cell growth and survival depend on a delicate balance between energy production and synthesis of metabolites. Here, we provide evidence that an alternative mitochondrial complex II (CII) assembly, designated as CIIlow, serves as a checkpoint for metabolite biosynthesis under bioenergetic stress, with cells suppressing their energy utilization by modulating DNA synthesis and cell cycle progression. Depletion of CIIlow leads to an imbalance in energy utilization and metabolite synthesis, as evidenced by recovery of the de novo pyrimidine pathway and unlocking cell cycle arrest from the S-phase. In vitro experiments are further corroborated by analysis of paraganglioma tissues from patients with sporadic, SDHA and SDHB mutations. These findings suggest that CIIlow is a core complex inside mitochondria that provides homeostatic control of cellular metabolism depending on the availability of energy.

Adverse effects of Hif1a mutation and maternal diabetes on the offspring heart.

BACKGROUND: Epidemiological studies show that maternal diabetes predisposes offspring to cardiovascular and metabolic disorders. However, the precise mechanisms for the underlying penetrance and disease predisposition remain poorly understood. We examined whether hypoxia-inducible factor 1 alpha, in combination with exposure to a diabetic intrauterine environment, influences the function and molecular structure of the adult offspring heart. METHODS AND RESULTS: In a mouse model, we demonstrated that haploinsufficient (Hif1a+/-) offspring from a diabetic pregnancy developed left ventricle dysfunction at 12 weeks of age, as manifested by decreased fractional shortening and structural remodeling of the myocardium. Transcriptional profiling by RNA-seq revealed significant transcriptome changes in the left ventricle of diabetes-exposed Hif1a+/- offspring associated with development, metabolism, apoptosis, and blood vessel physiology. In contrast, both wild type and Hif1a+/- offspring from diabetic pregnancies showed changes in immune system processes and inflammatory responses. Immunohistochemical analyses demonstrated that the combination of haploinsufficiency of Hif1a and exposure to maternal diabetes resulted in impaired macrophage infiltration, increased levels of advanced glycation end products, and changes in vascular homeostasis in the adult offspring heart. CONCLUSIONS: Together our findings provide evidence that a global reduction in Hif1a gene dosage increases predisposition of the offspring exposed to maternal diabetes to cardiac dysfunction, and also underscore Hif1a as a critical factor in the fetal programming of adult cardiovascular disease.

Asymmetric distribution of biomolecules of maternal origin in the Xenopus laevis egg and their impact on the developmental plan.

Asymmetric cell division is a ubiquitous feature during the development of higher organisms. Asymmetry is achieved by differential localization or activities of biological molecules such as proteins, and coding and non-coding RNAs. Here, we present subcellular transcriptomic and proteomic analyses along the animal-vegetal axis of Xenopus laevis eggs. More than 98% of the maternal mRNAs could be categorized into four localization profile groups: animal, vegetal, extremely vegetal, and a newly described group of mRNAs that we call extremely animal, which are mRNAs enriched in the animal cortex region. 3'UTRs of localized mRNAs were analyzed for localization motifs. Several putative motifs were discovered for vegetal and extremely vegetal mRNAs, while no distinct conserved motifs for the extremely animal mRNAs were identified, suggesting different localization mechanisms. Asymmetric profiles were also found for proteins, with correlation to those of corresponding mRNAs. Based on unexpected observation of the profiles of the homoeologous genes exd2 we propose a possible mechanism of genetic evolution.

The role of nitric oxide during embryonic epidermis development of Xenopus laevis.

Nitric oxide (NO) is a potent radical molecule that participates in various biological processes such as vasodilation, cell proliferation, immune response and neurotransmission. NO mainly activates soluble guanylate cyclase, leading to cGMP production and activation of protein kinase G and its downstream targets. Here we report the essential role of NO during embryonic epidermis development. Xenopus embryonic epidermis has become a useful model reflecting human epithelial tissue composition. The developing epidermis of Xenopus laevis is formed from specialized ionocytes, multi-ciliated, goblet and small secretory cells. We found that NO is mainly produced in multi-ciliated cells and ionocytes. Production of NO during early developmental stages is required for formation of multi-ciliated cells, ionocytes and small secretory cells by regulation of epidermal-specific gene expression. The data from this research indicate a novel role of NO during development, which supports recent findings of NO production in human mucociliary and epithelium development.